Characterisations of Classical and Non-classical states of Quantised Radiation 
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A new operator based condition for distinguishing classical from non-classical states of quantised 
radiation is developed. It exploits the fact that the normal ordering rule of correspondence to go from 
classical to quantum dynamical variables does not in general maintain positivity. It is shown that 
the approach naturally leads to distinguishing several layers of increasing nonclassicality, with more 
layers as the number of modes increases. A generalisation of the notion of subpoissonian statistics 
for two-mode radiation fields is achieved by analysing completely all correlations and fluctuations 
in quadratic combinations of mode annihilation and creation operators conserving the total photon 
number. This generalisation is nontrivial and intrinsically two-mode as it goes beyond all possible 
single mode projections of the two-mode field. The nonclassicality of pair coherent states, squeezed 
vacuum and squeezed thermal states is analysed and contrasted with one another, comparing the 
generalised subpoissonian statistics with extant signatures of nonclassical behaviour. 

I. INTRODUCTION 

Electromagnetic radiation is intrinsically quantum mechanical in nature. Nevertheless it has been found extremely 



Oh: 

' , . 

fruitful, at both conceptual and practical levels, to designate certain states of quantised radiation as being essentially 
^ "classical" , and others as being "non-classical" It is the latter that show the specific quantum features of ra- 
^ diation most sharply. Some of the well known signs of nonclassicality in this context are quadrature squeezing 
^^antibunching and subpoissonian photon statistics Q . 

ILJ The purposes of this paper are to present a new physically equivalent way of distinguishing classical from nonclassical 
. ^ States of radiation, dual to the customary definition and based on operator properties; to point out the existence of 
several levels of classical behaviour, with a structure that gets progressively more elaborate as the number of modes 

?-j increases; and finally to give a complete discussion of signatures of nonclassical photon statistics for two-mode fields, 
. . . working at the level of fluctuations in photon numbers. 

The contents of this paper are arranged as follows In section II we develop a criterion based on operator expectation 
values, to distinguish between classical and nonclassical states of radiation. The basic idea is that the normal ordering 
rule of correspondence between classical dynamical variables and quantum operators, while being linear and translating 
reality into hermiticity, does not respect positivity. If this potential nonpositivity does not show up in the expectation 
values of operators in a certain state, then that state is classical; otherwise it is nonclassical. Section III explores this 
new approach further and shows that, as the number of independent modes increases, the classification of quantum 
states gets progressively finer; several levels of nonclassicality emerge. This is shown in detail for one and two mode 
fields, and then the trend becomes clear. Section IV analyses in complete detail the properties of two mode photon 
number fiuctuations, stressing the freedom to choose any normalised linear combination of the originally given modes 
as a variable single mode. The well known Mandel parameter criterion for sub-Poissonian statistics for a single mode 
field is extended in full generality to a matrix inequality in the two-mode case. It is shown that certain consequences of 
this inequality transcend the set of all single mode projections of it, and are thus intrinsically two-mode in character. 
Explicit physically interesting examples of this situation are provided, and the well known pair coherent states are 
also examined from this point of view. Section V presents some concluding remarks. 
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II. THE DISTINCTION BETWEEN CLASSICAL AND NON-CLASSICAL STATES - AN OPERATOR 

CRITERION 



We deal for simplicity with states of a single mode radiation field, though our arguments generalise immediately to 
any number of modes. The photon creation and annihilation operators and a obey the customary commutation 
relation 



\a, cl'] = a — a ^ 1 ■ 



(2.1) 



The coherent states \z) are right eigenstates of d with (a generally complex) eigenvalue z; they are related to the 
states |n) of definite photon number (eigenstates of a) in the standard way: 
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(2.2) 

A general (pure or mixed) state of the one-mode field is described by a corresponding normalised density matrix p: 

^ p>0, Tr p = 1 • (2.3) 
It can be expanded in the so-called diagonal coherent state representation Q : 



— 4>{z) \z){z\, 

IT 



While hermiticity of p corresponds to reality of the weight function ( 
quantity, namely a distribution of a well-defined class. 



(2.4) 

z), the latter is in general a singular mathematical 

Namely, p is 



The conventional designation of p as being classical or nonclassical is based on the properties of M^) 
{z) is everywhere non- negative and not more singular than a delta function ||]: 



said to be classical if i 



p classical 
p nonclassical 4^ 



(z) > 0, no worse than delta function 



(2.5) 



It is clear that the conditions to be classical involve an infinite number of independent inequalities, since 0(z) > 
has to be obeyed at each point z in the complex plane. This is true despite the fact that the condition p > means 
that the "values" of (j>{z) at different points z are not quite "independent" . One realises this by recalling that every 
classical probability distribution over the complex plane is certainly a possible choice for 0(z), with the corresponding 
p being classical. 

The above familiar definition of classical states deals directly with p and (t>{z). Now we develop a d ual, but 
equivalent, definition based on operators and their expectation values. As is well known, the representation (2.4) for 
p is closely allied to the normal ordering rule for passing from classical c-number dynamical variables to quantum 
operators. Within quantum mechanics we know that an operator F is completely and uniquely determined by its 
diagonal coherent state matrix elements (expectation values) {z\F\z). Moreover, hermiticity of F and reality of (z|F|z) 
are precisely equivalent. Any (real) classical function f{z*,z) determines uniquely, by the normal ordering rule of 
placing d^ always to the left of d after substituting z ^ d and z* ^ d^, a corresponding (hermitian) operator F^ as 
follows: 

Normal ordering rule 
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{z\Fn\z) = f{z\z), 

f real <^ Fat hermitian (2-6) 



The connection with the representation (2.4) for p is given by 

Tr 



{pFn) = J ^<t>{z) f{z\z)- (2.7) 
It is an important property of the normal ordering rule that, while it translates classical reality to quantum hermiticity, 



it does not preserve positive semidefiniteness. More explicitly, while by eq. (2_^) -F/v > implies f{z*,z) > 0, the 
converse is not true. Here are some simple examples of nonnegative classical real f{z*,z) leading to indefinite 
hermitian Ff^: 

, 2 



f{z*,z) = {z* + zf ^FN={a^ + a) -1; 
f{z*,z) - (z* + zf ~^Fn= ({o) + a)' - 3)' - 6 ; 

/(z*,z) = e-^*^ ^ ^z*^z^^Fn = Y. Cn\n){n\- (2.8) 

n=0 ' ri=0 

In the last example, the real constants C„ can certainly be chosen so that some of them are negative while maintaining 
f{z*i z) > 0; this results in Fn being indefinite. 

We thus see that when the normal ordering rule is used, every _Fjv > arises from a unique f{z*, z) > 0, but some 
(real) f(z*,z) > lead to (hermitian) indefinite F]\[. So in a given quantum state p, the operator Fn corresponding 
to a nonnegative real classical /(z*, z) could well have a negative expectation value. // this never happens, then p is 



classical. That is, as we see upon combining eqns.( 2^) (2/7): if for every /(z*,z) > the corresponding Fn has a 



nonnegative expectation value even though _Fjv may be indefinite, then p is classical. Conversely, p is nonclassical if 
there is at least one f{z*, z) > which leads to an indefinite Fn whose expectation value is negative. 

We can convey the content of this dual operator way of defining classical states also as follows: while the normal or- 
dering rule allows for the appearance of "negativity" in an operator Fjv even when none is present in the corresponding 
classical f{z*, z), in a classical state such negativity never shows up in expectation values. 

Purely by way of contrast, we compare the above with what obtains when the antinormal ordering rule - substituting 
z a, z* — > followed by placing a to the left, to the right - is used to pass from classical f{z*, z) to quantum 



F ||7|1. In place of eqns.( |2J) (^ we have: 
Antinormal ordering rule 



. . - [ d z 

real/(z*,z) —> hermitian Fa = / f{z*,z) \z)(z\ , 

J TT 

[pFa) = / ^ {z\p\z) J{z\z)- (2.9) 



Tr 

Now the situation is that f{z*,z) > certainly implies Fa > 0, but some Fa > arise from indefinite classical 
f(z*,z). A simple example is: 

f{z*,z) = (z* + z)2 - I ^ Fa = (at + a)' • (2.10) 

Thus the function and operator relationships are opposite to what we found in the normal ordering case. The direct 
characterization of classical states (whose definition in any case is a convention based on the normal ordering rule) is 
however not very convenient with the antinormal ordering convention. 

In passing we may note an interesting aspect of the Weyl or symmetric rule of ordering 1^, and the associated 
Wigner distribution description of a state p 9| . Here it is more convenient to deal with the real and imaginary parts 
of z and a: 

z = ^{q + ip), z*^^{q~ip); 

a = -^{q + ip), a^ = ^{q-ip), 

[9, P] = z • (2.11) 
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The Weyl ordering rule maps single classical exponentials into single operator exponentials, 



exp(i X q + i ^ p) ^ exp(i X q + i ^ p) , X and real ; (2-12) 

and then extends this by linearity and Fourier transformation to general functions: 
Weyl rule 



C30 oo 

f{q,p) = J dX J d^i f{X,^i)exp{i Xq + i fip) 



-OO — OO 

OC OO 



Fw — dX dfi f{X, fi) exp{i X q -\- i fi p) , 



-OO — OO 



/ real Fw hermitian (2.13) 
The Wigner distribution W{q^p) for a state p is given in terms of the configuration space matrix elements of p: 

OO 

W{q,p) = / dq'iq - ^q'\p\q+lq') ^ (2.14) 



and the rule for expectation values ties together eqns. (2.13) (2.14) 



OO OO 



Tl-i^pFw)^ J J dqdpW{q,p) f{q,p)- (2.15) 

— OO — OO 

From the algebraic point of view the Weyl rule stands "midway" or symmetrically between the normal and the 
antinormal ordering rules. This may lead one to hope that it removes the mismatch and resolves the problem of 
preserving positivity in both directions in passing between f{q,p) and Fw- However this does not happen at all. As 
is well known there are nonnegative f{q,p) yielding indefinite Fw, and nonnegative Fw leading back to indefinite 
f{q,p)- Here are simple examples: 

f{q,p) = 6(q)6{p) — > Fw = parity operator with eigenvalues ± 1 ; 

Fw = |1)(1| ^ f{q,p) = -(q'+p'- I:) cxp ( V - p') ' (2.16) 



TT V 2, 

While these remarks illuminate in terms of operator properties the relations among the three ordering rules, the 
classification of states p into classical and nonclassical ones is based most simply on the normal ordering rule. It is 
clear that all these considerations extend easily to any number of modes of radiation. 

While our operator based approach to the identification of nonclassical states is conceptually complete, we note 
that most of the extant criteria of nonclassicality involve a simple extension of our formalism. Namely, one often 
has to consider nonlinear functions of expectation values of several operators, which cannot be simply expressed as 
the expectation value of a single state-independent or autonomous operator. This is always so when one deals with 
fluctuations. Thus quadrature squeezed non-classical states are defined via inequalities involving the fluctuations 
(Ag)^ and (Ap)^. Similarly amplitude squeezing or sub-Poissonian statistics deals with the fluctuation in photon 
number. These remarks apply also to other criteria of non-classicality such as higher order squeezing [ p^ and the one 
based on matrices constructed out of factorial moments of the photon number distribution The one interesting 
exception to these remarks and which is fully covered by our formalism, is the case of antibunching . Here one is 
concerned with the expectation values of the single time-dependent operator 

F{t, t + T)=: i{t + T)i{t) : - : /(t)/(t) : (2.17) 

where I{t) is the intensity operator at time t. In this sense, the criterion for antibunching is qualitatively different 
from the other familiar ones. 
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III. LEVELS OF CLASSICALITY 



A. The single mode case 

We begin again with the single mode situation and hereafter deal exclusively with the normal ordering prescription. 
(Therefore the subscript N on Fjv will be omitted). Suppose we limit ourselves to classical functions f{z*,z) which 
are real, nonnegative and phase invariant, that is, invariant under z e*"z. An independent and complete set of 
these can be taken to be 

U{z*, z) = e-^*^ z*^z^/n\ , n = 0, 1, 2, . . . , (3.1) 
since they map conveniently to the number state projection operators: 

/„(z*,z)^f(") = , 71 = 0,1,2,... (3.2) 

A general real linear combination f{z*,z) = ^ C„ /„(z*,z), even if nonnegative, may lead to an indefinite F, as 



seen at eqn. (p.Sj). 

If we are interested only in the expectation values of such variables, we are concerned only with the probabilities 
p{n) for finding various numbers of photons; for this purpose an angular average of (j){z) is all that is required: 



p{n) = {n\p\n) = Tr ip\n)(n\) 

(z) e-"'"z™^"/n! 



TT 



dl P{I) e-' /"/n! , 



Now while (/)(z) > certainly implies P{I) > 0, the converse is not true. Thus one is led to a three-fold classification 
of quantum states p p^ : 

p classical ■^=> (f>{z) > , hence P{I) > ; 
p semiclassical P{I) > but (j){z) ^ ; 
p strongly nonclassical <;=^ P{I) ^ , so necessarily (/)(z) ^ • (3-4) 



The previous definition (2.5) of nonclassical p based on 0(z) alone is now refined to yield two subsets of states, the 
semiclassical and the strongly nonclassical. The semiclassical states do have the following property: 

oo 

p semiclassical =^ Tr (^pF^ > if f{z*, z) = ^ C„ /„(z*, z) > • (3.5) 

ri=0 

However, in addition, there would definitely be some phase noninvariant f{z*,z) > for which F is indefinite and 
Tr (/5f) < 0. It is just that this extent of nonclassicality in p is not revealed by the expectation values of phase 

invariant variables, or at the level of the probabilities p{n) [p^ . 

It is clear that the classification (3.4) is U{1) or phase invariant. That is, p retains its classical, semiclassical or 
strongly non-classical character under the transformation (/>(z) — s- 4>'{z) = (j}{ze'^°') 

As examples of interesting inequalities obeyed if p is either classical or semiclassical, we may quote the following 
involving the factorial moments of the photon number probabilities p{n) : 

7™ = Tr(pat'"a™) 

— 4>{z) {z*zY' 

TT 
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= J d/p(/)/" 



= ^ n!/(n-TO)! > 0, m = 0, 1, 2, . . . ; 

n—m 

p classical or semiclassical P{I) > =^ 7™ 7„ < 7„i+„ < V72m72n (3.6) 

Violation of any one of these ine qual ities implies p is strongly nonclassical. 

The inequalities quoted in eq. (3^) above clearly involve an infinite subset of the photon number probabilities p{n). 
However one can easily construct far simpler inequalities involving a small number of the p(n)'s, violation of any of 
which also implies that p is strongly non-classical. For example, from eqs. (3.1) (3.2), for any nonnegative integer uq 
and any real a, 6 we have the correspondence 

f{z\ z) = e-"*" ^ (a + bzz*f ^ 

tiq! 

F = a^\no){no\ + 2{no + l)ab\no + l){no + 1| + (no + l)(no + 2)b^\no + 2){no + 2| (3.7) 
Here f{z*, z) is nonnegative while F is indefinite if ah < 0. We then have the result : 

p classical or semiclassical P(/) > ^ 
a^p(no) + 2(^0 + l)abp{no + 1) + (^o + l)(no + 2)b^p{7io + 2) 



= ^ f dIP(I)e-^r°{a + blf > (3.8) 
no! J 



So again, violation of any of these " local" inequalities in p{n) implies that p is strongly nonclassical. 

A physically illuminating example of the distinction between classical and semiclassical p, and passage from one 
to the other, is provided by the case of the Kerr medium. The argument is intricate and rests on two well known 
results. The first is Hudson's theorem : if a (purestate) wavefunction V'o(9) has a nonnegative Wigner function 
Wo{q,p), then ipo{q) is Gaussian and conversely; in that case Wo((j',p) is also Gaussian. The second result is the 
general connection between ^(-z) and W{q,p) for any p: 



W{q,p) = 2 I ^e-2|---'l>(z'), z = —{q + ^p)- (3.9) 
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This means that for classical p with 0(z) > 0, W{q,p) > as well. Now imagine a single mode radiation field in an 
initial coherent state \zo) with 0o(^) = nS^^'^lz — zq), incident upon a Kerr medium [p^ . This initial state is pure, 
classical, and has a Gaussian wave function ipaiq)- The Kerr medium Hamiltonian is of the form 

-ffj^gj-j. = a a^a + /3(a^a)^ (3.10) 

Clearly the number states |n) are eigenstates of this Hamiltonian. Therefore the Poissonian photon number distribu- 
tion 

p{n) = e-'»IS/n\, h = z*z^, (3.11) 

of the input state jzo) is preserved under passage through the Kerr medium. Likewise the funtion P{I) = S{I—Io) is left 
unaltered. Therefore the output state lip), which of course is pure, is either classical or semiclassical. However the form 

of ^^j^err shows that the output wavefunction is non Gaussian^ Therefore by Hudson's theorem the corresponding 

M^(g,p) must become negative somewhere. Therefore by eq. ( |3.9| ) the output 0(z) cannot be nonnegative. Thus 
passage through the Kerr medium converts an incident coherent state, which is classical, into a final state which is 
semiclassical. 



B. The Two- mode case 



Now we sketch the extension of these ideas to the two-mode case. Here the operator commutation relations, number 
and coherent states, and the diagonal representation for p, are as follows: 
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[or, al] = Srs, [cLr, As] = [4' «1] = 0' '',5 = 1,2; 
|ni,n2) = (niln^!)-'/' (al)"' (4)"' |0, 0) , 

(^a\ai or 0^02^ \ni,n2) ~ (ni or 7^2) 712) ; 

k) = ki,Z2>=exp --z"fz ^ — rT72l"i''^2) 



J d^U) 4>iz)\z){z\ , 

d^i{z) = cPzi Sz^jif-- (3.12) 

It is convenient at this point to go when necessary beyond purely real classical functions /(z) in applying the normal 
ordering rule to obtain corresponding operators. From the general number states matrix elements of p we read off 
some operator correspondences generalising eqns. (3.1) (3.2): 

(n3,n4|/5|ni,n2) = Tr(/5|ni, 712) (^3, 714!) 

dfi{z) (l){z) e' 



ni!7T,2!7T,3!7i4! 

. *'ni *n2 n4, 

e-^-^- \ ^ ^ ^ ™> |nl,7^2)(r.3,n4| • (3.13) 

For one mode the phase transformations form the group UiX). For two modes this generalises to the group C/(2) of 
(passive) transformations mixing the two orthonormal single photon modes. At the operator level this means that the 
annihilation operators a,, experience a general ?7(2) matrix transformation conser ving total photon number. Moreover 
these transformations are unitarily implemented on the two-mode Hilbert space |16| : 

u = (Tirs) e C^(2) : 

lA(u) hr U(u)^^ = Usr , 

Uiu) a\ lA{u)^^ — u*j. al , 
W(M)W(ti)t= 1, 

U{u')U{u)= U{u'u)- (3.14) 

For later use we give here the actions of these unitary operators hl{u) on monomials formed out of al and a^, on the 
number states and on coherent states: 



'a e C/(2), a e SU{2) : 
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\/ U + my.(j - m)l ' ' ' + m')l{j -rn')l 

'•j-\-'m " j — m - j-\-rn ■^j — rn 

V(j +m)!(j - m)! ^/{j +m'y.{j -m'y. 

U{u)\j +m,j~m)= ^ Dl^]ja)\j +m',j- m') , 

771 ' 

j^O, 1/2, 1,..., 
m, m' ^j, j - 1, . . . , -j ; 
U{u)\z) = \u* z) ■ (3.15) 

We have chosen the exponents of a's and a^'s, the numerical factors, and the number operator eigenvalues, in such a 
way that the results can be expressed neatly using the SU (2) representation matrices in various unitary irreducible 
SU{2) representations, namely the 2?-functions of quantum angular momentum theory [TtI] . 

To motivate the existence of s everal layers of classicality, we now generalise the single mode [/(l)-invariant real 
factorial moments 7„ of eqn. (3.6) to two-mode quantities which conserve total photon number and also transform in 
a closed and covariant manner under SU{2). For this purpose, keeping in mind eqns. ( 3.15| ), it is convenient to start 
with the (in general complex) classical monomials 
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/mim2(— ) — ^jmim^^l 



*j+mi ^j—mi j+m2 



^2 1 



Nj,n,,n., = [{i + mi)! (j - mi)! (j + mz)! {j - ma)! t^'"^ , 

] = 0, 1/2, 1, . . . , 
mi, 1712= j, • 



(3.16) 



The total power of z is equal to that of z} , hence these are ?7(l)-invariant. The correspondmg operators and then 
SU (2) transformation laws are: 



Jmim2 \— '— / ^mim2 ~ -'''jmim2"l "2 "l "2 i 

aeS'C/(2) :ZY(a)F;;, „, Zi(a)-^ = V D^^] {a)D^'] {a)*F\ ,■ 



(3.17) 



For a given two-mode state p we now generalise the factorial moments 7„ of eqn. (3.6) to the following three-index 
quantities: 



^(J) = Tr f o 



N- Tr 



J-ni2 
2 



(3.18) 



Their S*?/ (2) transformation law is clearly 



p =U{a) pU{a)-^ : 

m'^ ,7712 



(3.19) 



In the last line for each fixed j the generalised moments 7rni7ra2 have been regarded as a (hcrmitian) matrix of dimension 
(2j + 1). 

On account of the fact that the total photon number is conserved in the definition of these moments, calculation of 
7mi77i2 does not require complete knowledge of 0(z) but only of a partly angle averaged quantity ^(/i, /2, &)'■ 



2ir 



de 

2^ 



■p(/i,/2,6')(/i/2)^(/i//2)^/^'"'+"'^ • e*^™!-"^)'^ 



2tt 



7'(/i, /2, ^) = / |i <^ (/y V«M2^/V^^+^)) 



(3.20) 



It is clear that these moments 7ml7Ti2 involve more than just the photon number probabilities p{ni-, ^2) which are just 
the "diagonal" case of the general matrix element in eqn. (3.13): 



p{ni,n2) = {ni,n2\p\ni,n2) 

00 00 

= Jdh j dl2P{h,l2)e-''-'-r^'ri'/n^\n2\ 
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P{h.l2)= j ^V{Iul2,t 



(3.21) 



This is the two-mode version of eqn. (3.3). The subset of "diagonal" moments 7mm are calculable in terms of p(ni, 77.2) 
or P(/i,/2): 
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oo oo 



= jdh J dh P(/i,/2)/f+"/r™/(j + {j my 



= p{ni,n2) nil n2l/ {fii — j — my.{n2 ~ j + my.{j + my.{j my. (3.22) 

ni ,712 

However under a general SU{2) mixing of the modes, the expressions 7mm, pini,n2), P(/i,/2) do not transform 
in any neat way among themselves, and one is obliged to enlarge the set to include the more general 7mlm2 and 
7'(/i, J2,0). (In particular, for these, the probabilities p{ni,n2) are inadequate). When this is done we see the need 
to deal with both the quantities V{Ii,l2, P{Ii,l2) derived from (j>{z) by a single or a double angular average. One 
can therefore distinguish four levels of classicality for two-mode states: 

p classical <^ (/)(z) > 0, (hence ^(/i, /a, 6*), P(/i, ^2) > 0) ; 
p semiclassical / <^ P(/i, /2, 6*) > (hence P(/i, 12) > 0) , but <p{z) ^ ; 
p semiclassical // <^ P(/i, 12) > but V{Iij2,d) ^0 (hence (/)(z) ^ 0) ; 
p strongly nonclassical ^ P(/i, 12) ^ 0(hence (/)(z), V{h,l2,0) ^ 0)- (3.23) 

These definitions can be cast in dual operator forms. For example, for semi classical -I states, we can say that for 
any classical real nonnegative overall U{1) phase invariant / {z^ ,z}j the corresponding operator F has a nonnegative 
expectation value, while this fails for some / {z} outside this class. In the semi-classical-II case, we have to further 
limit / i^z^ , to be real nonnegative and invariant under independent U{1) x U{1) phase transformations in the two 

modes, to be sure that the expectation value of F is nonnegative. 

At this point we can see that these levels of classicality possess different covariance groups. Since under a general 
U{2) transformation Zi(u), M g U{2), the function 0(2) undergoes a point transformation, 0(z) — > (l)'{z) = (f>{u^z_), we 
see that the property of being classical is preserved by all U{2) transformations. On the other hand, a general U{2) 
transformation can cause transitions among the other three levels. The point transformation property is obtained 
for V{Ii,l2,0) and P(/i,J2) only under the diagonal U{1) x U{1) subgroup of U{2); in fact P(/i,/2) is invariant 
under U{1) x U{1), while V{Ii, l2,0) suffers a shift in the angle argument 9. Thus one can see that each of the three 
properties of being semiclassical-I, semiclassical-II or strongly non-classical is only U{1) x J7(l) invariant. 

As the number of modes increases further, clearly the hierarchy of levels of classicality also increases. 

Generalising inequalities of the form (3.6) for the diagonal quantities -jmln, for p classical or semiclassical-I or 



semiclassical-II, is quite straightforward, since then we deal with the two modes separately. The more interesting. 



and quite nontrivial, problem is to look for matrix generalisations of eqn. (3.6), bringing in the entire matrices 
7^^^ — (^^m{m2^, and looking for inequalities valid for states of the classical or semi-classical-I types. (Of course for 

any quantum state p we have the obvious property that 7^-''', for each j, is hermitian positive semidefinite. This is the 
two-mode generalisation of 7„ > in the one- mode case) . However this is expected to involve use of the Racah-Wigner 
calculus for coupling of tensor operators, familiar from quantum angular momentum theory, inequalities for reduced 
matrix elements, etc Jl8{ . 

In the next Section we undertake a study of the particular case j = 1 which involves at most quartic expressions in 
a's and d^'s. This is just what is involved in giving a complete account of the two- mode generalisation of the Mandel 
Q-parameter familiar in the single-mode case. 

IV. GENERALISED PHOTON-NUMBER FLUCTUATION MATRIX FOR TWO-MODE FIELDS 

For the one-mode case, with the single photon number operator JV — d^a, we have some ob viou s inequalities valid 
in all quantum states, and others valid in semiclassical and classical states as defined in eqn. ( ^.4[ ): 

Any state : 

(iV) = Tr (piV) = 71 > ; (4.1a) 

( : iV^ :) = Tr (pa^'^a^) = 72 > ; (4.1b) 

{N^) = Tr [pa^aa^a) = 72 + 71 > ; (4.1c) 

(A7V)2 = (iV2) - (iV)2 = ((iV - (iV))2) = 72 + 71 - 7? > ; (4.1d) 



9 



semiclassical or classical state 

(: :) - (7V)2 = (ANf - {N) = 72 - 7? > • (4.1e) 

(Here the dots : : denote normal ordering). The Mandel Q-parameter is defined as |^ 

^^(AiVf-W^T^^ (4.2) 
{N) 71 

and it has the property of being nonnegative in classical and semiclassical states. Conversely if Q is negative, the state 
is definitely strongly nonclassical. The two cases Q > and Q < correspond respectively to super and subpoissonian 
photon number distributions. 

The inequalities (4.1) are not all independent, as some imply others. We now give the generalisation of these in 
matrix form, to two-mode states. 

We have to deal with four independent number-like operators A^^, /i = 0, 1, 2, 3 which we define thus: 

alas = ^(crp)rs^^r (4.3) 

(Here ctq and aj are the unit and the Pauli matrices, and the sum on fi goes from to 3). The expectation values of 
Nfj_ in a general state p are written as : 



alas) 



1 



Tr (palas) = ■^icrii)rs n^,. (4.4) 



Thus and the matrix 7(^/2) = {^mim^ are essentially the same. Since the 2x2 matrix ((aja^)) is always hermitian 
positive semidefinite, we see that the generalisation of inequality ( 4.1a ) to the two-mode case is 

no - Izil > • (4.5) 

(All components of are real). It may be helpful to remark that the matrix 7^^/^) is analogous to the coherency 
matrix, and the quantities to the Stokes parameters, in polarisation optics |19| . 

Now we consider quadratic expressions in which are upto quartic in and combined. To handle their normal 
ordering compactly, we first define certain quadratic expressions in ar, and their hermitian conjugates: 

Aj = i a^ a2 Oja , Al^ = —i a'(jj(T2a* , j = 1, 2, 3 ; 

ar&s = icrj<T2)^s ' ^l^l = ^i'^2Crj)rs • (4.6) 

Under the action of the unitary operators (a) representing SU{2), both Aj and Aj transform as real three-dimensional 

Cartesian vectors. Now we can easily express the result of writing the product N^Ni, as a leading normally ordered 
quartic term plus a remainder: 

iV^TV^ = : TV^iV^ : 4- {(.^^x + i eop,.A) ^Va , 
: N^,N^ : = t^,^jkA]Ak , 

<^tiux = Sf^iySxo + Sf,oS,yx + St^aSfj^x - 2 J^o (^i/o Sxo- (4.7) 

Here ea-^ux is the four-index Levi-Civita symbol with 60123 = 1- So the anti-commutators and commutators among 
Nf^ and Ni, are: 

i {N^, N,} = t^,,jkA]Ak + (if^uxNx , (4.8a) 
N,] = 2 I eo^uxNx (4.8b) 
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(These latter are just the C/(2) Lie algebra relations). To accompany n^, in a general state we denote the expectation 
values of ^]^fc by Qjk- 



Clearly, (qjk) is basically the matrix 7 



{A]Ak 



Tr 



{pA]Ak) 



qjk ■ 



(4.9) 



statement is the generalisation of inequality ( 4.1b| ). We can also generalise the inequalities (4.1c) (4. Id) by saying 
that for any quantum state the two matrices with elements given by 

(i{iV^, N^}) = i^^jTc^jfe + ^M^A^A , (4.10a) 

A{%,N,) EE i({7V^ - {%) , - (iV,)}) 

= tf^^jkqjk + ifi,y\n\ - (4.10b) 

are both 4x4 real symmetric positive semidefinite. As in the one-mode case the inequality obeyed by ^A^^, ^i')) 

implies the one obeyed by the anticommutator matrix ^(^ |-/Vp, 7Vy|)^ 

Now we search for matrix inequalities which are valid in two-mode classical or semi-classical-I states, but not 
necessarily in semi-classical-II or strongly nonclassical states. The key ingredient is the formula 



and is always a 3 x 3 hermitian positive semidefinite matrix. This 



1 



{N^){N.) 



00 oc 27r 



jl/2 

jl/2 ^.0 



(4.11) 



We can now draw the following conclusion: 
Classical or Semi-classical-I state 



1 



= (A(7V^,iV,) 



> 0' 



(4.12) 



This is the intrinsic two-mode expression of super-poissonian statistics, and its violation (possible only in semi-classical- 
II or strongly nonclassical states) is an intrinsic signature of two-mode subpoissonian photon statistics. What makes 
this criterion nontrivial is the fact that for any obeying eq. ( [1.5| ) the 4x4 matrix {If^^xux) is real symmetric 
positive semi-definite. 

It is interesting to pin down the way in which this matrix inequality ( 4.12 ) can go beyond a single- mode condi- 
tion pc[ | . The most general normalised linear combination of the two mode-operators is determined by a complex 
two-component unit vector a: 

a{a) = a}a = a*dr , 
a' a = 1 ; 
[a{a) , a(a)^] = 1 • 



(4.13) 



For every such choice of a single mode, the inequality ( [4.12| ) does imply the single- mode inequality (l.le). We can 
see this quite simply as follows. Given a, we define the real four-component quantity ^^(a) by 



Co (a) = le(a)| = 1/2- 



(4.14) 
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Then, using the completeness of cr^i expressed by 

(o-p)rs i<ytj.)tu = 2 Sru Sgt , (4-15) 

we have the consequences: 

£.^^{a)N^. = a(a)ta(a) = N{a) , 

e^..xU^DU(£) = U^i)- (4.16) 

Indeed we easily verify that (leaving aside = identically) 

If^ux^t^^u = 6 =^ either = 1^1 = 1/2 



1 

some a obeying a' a = 1 , 

Co = l,e = 0- (4.17) 



Saturating the left hand side of ( 4.12 ) with the latter possibility, = (5^0, leads to the superpoissonian condition for 
the total photon number distribution. Saturating it with (a) we get as a consequence: 

(AiV(a))2 - (iV(a)) > 0, any a- (4.18) 

In this way the two-mode matrix "superpoissonian" condition ( [4.12 ) implies the scalar single mode super poissonian 
condition (4.1e) for every choice of normalised single mode with annihilation operator a{a), as well as for the total 

photon number. 

However, it is easy to see that the infor mation contained in the matrix inequality (4-1^) *s not exhausted by the 
collection of single mode inequalities ( 4-.1S) for all possible choices of (normalized) a. Denoting the real symmetric 
matrix appearing on the lefthand side of ( 4.12| ) by (A^^), 

A^, = A{N^,N,)-£^,x{Nx) , (4.19) 

it is clear that 

Cm 6 > for aU obeying = 1^1 = 1/2 

¥> (A^u) > (4.20) 

Indeed, the lefthand side here reads in detail: 

C^A^^C^ = ^Aoo + Aoj^j + (.j£.kAjk ; (4.21) 

and the nonnegativity of this expression for all 3- vectors with |^| = 1/2 cannot exclude the possibility of t he 3 x 3 
matrix {Ajk) having some negative eigenvalues. Part of the information contained in the matrix condition (4.12) is 
thus irreducibly two-mode in character, a sample of this being: 

(A^.) >0^ {A,k) > • (4.22) 

Admittedly to a limited extent, this situation is analogous to some well known properties of Wigner distributions. 
Thus the marginal distributions in a single variable obtained by integrating W{q,p) with respect to p or with respect 
to q (or any real linear combination of q and p) are always nonnegative probability distributions, even though W{q,p) 
is in general indefinite. So also here, it can well happen that for a certain state both Aqq and ^fj.{a)Af^i,^^{a) are 
nonnegative for all a , yet (A^j,) is indefinite. 

There exists in the literature a well known inequality for two-mode fields, which when violated is a sign of nonclas- 
sicality [Q. It reads: 

(ni(ni - 1) + n2in2 - 1) - 2hih2) > 0, 

ni~a\ai, n2 — 02^2, (4.23) 

and evidently involves only diagonal elements of the matrix (A^^). After rearranging the operators in normal ordered 
form one can see that 
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(ni(ni - 1) + n2{n2 - 1) - 2nin2) 
— {a[ + — 2a[a2aia2) 



= 2(911 



A 



33 -1- 

{a\ai 



922 — 933 J 



(4.24) 



By our analysis, in any classical or semiclassical-I state the matix {A^^) is positive semidefinite, so in particular ^433 



and even more so the expression ^33 + rtg, are both nonnegative. Thus the inequality ( 4.23| ) is certainly a necessary 



condition for classical and semiclassical-I states. Conversely, if ( 4.23| ) is violated and ^33+^3 is negative, then certainly 



A^3 is negative as well, and the state is either semiclassical-II or strongly non-classical. However this condition is 
unnecessarily strong since it asks for A33 to be less than — n§; as we have shown, even the weaker condition A33 < 
is sufficient to imply that the state is semiclassical-II or strongly nonclassical. Vic e ver sa, our necessary condition 
^33 > for a classical state or semiclassical-I state is stronger than the condition (4.23). In both directions, then, 
our conditions are sharper than the ones existing in the literature. 



We conclude this Section by presenting a few examples bringing out the content of the matrix condition ( 4.12 ), in 
particular the possibility of its containing more information than all single-mode projections of it. 

(a) Pair-coherent states : 



These are simultaneous eigenstates of 0.102 and a\ai — 0)20,2 p^ ]: 

aia2|C,9) = C\C,,q), C e A 
{a\ai - 4a2)|C,9> = q\C,,q), q = 0,±1,±2- ■ 



(4.25) 



For q > these states are given by 

00 

\C,q)=N,J2 



C" 



n=0 



[n!(n-f g)!]i/2 



\n + q,n) 



(4.26) 



where Nq is a normalisation constant. It is known th at in these states the second mode already shows sub- 
poissonian statistics ||23|. Thus if we write the matrix ( 4.19 ) for these states as {Afj,i,{(,q)), then even without 
having to nontrivially mix the modes we find: 



1 



a= (0, 1)^ 



eM(a) = ^a'^fTMa=(l/2, 0, 0, -1/2 



(4.27) 



The matrix A{C,, q) is indefinite and the pair coherent states are therefore neither classical nor even semiclassical- 
I. Consistent with this, a direct numerical study of the least eigenvalue / {A{(, q)) of(^^j/(C, g)) for sample values 
of C and g, does show it to be negative. 

(b) Two-mode squeezed vacuum: It has been shown elsewhere [ p4| that a two mode squeezing transformation is 
characterised by two independent intrinsic squeeze parameters a and b obeying a > 6 > 0. A representative of 
such a transformation is 



Z^(o)(a,5) =exp 



(a-b) ^^2 .2^1 



exp 



(a + ^.~t2 .2^ 
(«2 -02) 



(4.28) 



The case a = b essentially corresponds to the second mode alone being squeezed. For general a ^ b we have 
genuine two-mode squeezing; while the (Caves-Shumaker) limit 6 = involves maximal entanglement of the two 
modes. We restrict our analysis to this limit in the sequel. Then the two-mode squeezed vacuum is characterised 
by the single parameter o and is 



Z^(°na,0)|0,0) = 



exp 



a /.t2 ,-2 



it2 



|0,0) 



(4.29) 



The matrix {A^i,{o)) can be explicitly computed and happens to be diagonal: 
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{Af,^{a)) = Diag ( ^ (-3 + 7 cosh(2 a)) Sinh(a)^ , 2 cosh(2 a) Sinh(a)^ , 

-2Sinh(a)^ , 2 cosh(2a) Sinh(a)^ 



(4.30) 



We see that for all a > this is indefinite, since the third eigenvalue A22{ a) is s trictly negative. This is displayed 
in Figure la, for a in the range < a < 1. Thus for all a > the state (4.29) is definitely neither classical nor 
semiclassical-I. On the other hand the leading diagonal element (eigenvalue) Aqq (a) dominates the others in the 
sense that for all choices of single mode the "expectation value" of A{a) is nonnegative: 



^tiia)Af,i,{a)^i,ia) >0, all a 



(4.31) 



Thus the squeezed vacuum (4.2E ) displays nonclassicality via subpoissonian statistics in an intrinsic or irreducible 
two-mode sense which never shows up at the one mode level for any choice of that mode. This is to be contrasted 
with the case of pair-coherent states discussed previously. At the same time the state ( 4.29| ) is also quadrature 
squeezed for all a > 0. Thus both these nonclassical features are present simultaneously. 

(c) Two-mode Squeezed Thermal state: This is defined as follows (we again limit ourselves to the case 6 = 0): 



poifi) = (1 - e-^f exp \-fi{a\ax + 



(4.32) 



At zero temperature /3 — > oo this goes over to the previous example (b). Once again the matrix (Ap,y(a, /?)) can 
be computed analytically and it turns out to be diagonal: 



{A,,,{a,l3))^{-l + e^f X 

I i(l3-14e'3-t-13e2/' + 20(l-e2'3) cosh(2 a) + 7 (1 + e'')^ cosh(4a)) , \ 
i((l- 6^3)2 + 2 (1-62/3) cosh(2a) -f (l-fe'5)2 cosh(4a)), 

l + e2^ + (l-e2'3) cosh(2a), 
i((l-e^')2 + 2(l-e2'3) cosh(2a) + (H-e'3)2 cosh(4a)) / 



Diag 



(4.33) 



Now the third element ^22(07/3) can become negative for low enough temperature T = or high enough 
squeeze parameter a. The variation of the least eigenvalue I (A(a, /?)) of A{a, /?) with respect to a in the range 
< a < 1, for various choices of /?, is shown in Figures (lb,c,d). One can see that if the temperature is not too 
high, for sufficiently large a the element A22(a,/3) becomes negative, indicating that the state has then become 
semi-classical-II or strongly nonclassical. (In comparison we recall that for quadrature squeezing to set in the 
parameter a must obey the inequality a > In Coth('|) ( ||2^) ) On the other hand as in example (b), the leading 
element Aoo(a, /?) again dominates the others in the sense that 

ep(a)A^.(a,/3)C.(a) > 0, alia (4.34) 



So once again, when ^22(0, /3) < 0, the subpoissonian statistics is irreducibly two-mode in character. In 
Figures (lb,c,d) we have also indicated the value of the squeeze parameter a at which quadrature squeezing sets 
in. It is interesting to see that, for the states described here, at each temperature, the irreducible two- mode 
subpoissonian statistics occurs before squeezing. Therefore ( limiting ourselves to low order moments of <^(^z)) 
there exists a range of squeeze parameter where the only visible nonclassicality is through such subpoissonian 
statistics. 

The more general squeezed thermal state 

/5(/3,a,6) =Zi("'(a,fe)Po(/3)Z^<°^(a,6)-\ (4.35) 



has qualitatively similar properties. Detailed numerical studies presented elsewhere [po| have shown that these 
states also do not show subpoissonian statistics at the one-mode level. On the other hand, direct search for the 
least eigenvalue of (A^,y(a, 5, /3)) reveals that, for suitable values of (i, a, b, this is negative. 

We thus have several instructive examples of the situation indicated by eq. (4.20) 
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V. CONCLUDING REMARKS 



We have presented a dual operator and expectation value based approach to the problem of distinguishing classical 
from non-classical states of quantised radiation, and thus brought out the significance of this classification in a new 
physically interesting manner. As the number of independent modes increases, this approach leads to finer and yet 
finer levels of nonclassical behaviour, in a steady progression. This has been followed up by a complete analysis of 
photon number fluctuations for two- mode fields, and a comprehensive concept of subpoissonian statistics for such 
fields going beyond what can be handled by techniques developed at the one-mode level. 

In a previous paper we have set up the formalism needed to examine the possibility of two-mode fields showing 
subpoissonian statistics at the one-mode level in an invariant manner, by following the variation of the Mandel Q- 
parameter as one continuously varies the combination of the two independent modes into a single mode. One can see 
through the work of the present paper that that preparatory analysis is a necessary prerequisite to be able to pinpoint 
the aspects of subpoissonian statistics which are irreducibly two-mode in character. Examples (b) and (c) at the end 
of Section IV brin g ou t this aspect vividly. 



The inequality (4.23) has been strengthened by our approach to a sharper criterion to distinguish various situations: 



Classical or Semiclassical-I A33 > 0; 
^33 < ^ Semiclassical-II or strongly nonclassical. (5-1) 



From eqns. (|4.19| , |4.24| ) we see that A33 has the following neat expression: 



A33 = (Ani)' - (ni) + (Ana)' - ("2) - 2A {m,n2) 

= {{fii - n2f) ~ {{hi - h2)f - {hi + n2> (5.2) 

It is thus expressible solely in terms of expectations and fluctuations of the orig inal (unmixed) mode number operators 
hi,h2 and their functions. One can now see easily, again from equation ( [f.lQ ), that the statements ( pTl| ) are part of 



a wider set of statements involving only expectations of functions of hi,h2'. 

Aoo = [aNoY -{No), 

Ao3 = A30 = A (7Vo,iV3) - {N3), 

A33 = (^AN^y ~{No), 

No = hi + h2 , N3 ^ hi - h2; (5.3a) 

Classical or semiclassical-I ( I > 0; 

\ ^30 ^33 / 

^00 ^03 j <; g Semiclassical-II or strongly nonclassical. (5.3b) 
^30 ^33 / 

All other inequalities involving matrix elements such as Aqi, ^02,^13 • ■ • involve "phase sensitive" quantities going 
beyond hi and 7^2 ■ 

Going back to the matrix A — {A^^), we see that from its properties we cannot immediately distinguish between 
the classical and semiclassical-I situations, or between the semiclassical-II and strongly non-classical situations. In 
both the former, A is positive semidefinite; while if A is indefinite, one of the latter two must occur. It would be 
interesting, for pair coherent states or squeezed thermal states for instance, to be able to see, when A is indefinite, 
whether we have a semi-classical-II or a strongly non-classical state, and whether this depends on and varies with the 
parameters in the state. 
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